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Within a wide class of ferromagnetic and antiferromagnetic systems, quantum 
tunneling of magnetization direction is spin-parity dependent: it vanishes for mag- 
netic particles with half-integer spin, but is allowed for integer spin. A coherent-state 
path integral calculation shows that this topological effect results from interference 
between tunneling paths. 
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The search for new systems which exhibit quantum phenomena at the mesoscopic scale 
has led to a great deal of activity in recent years on nanometer-size magnetic particles. In 
such particles it is possible for the electronic spins to be locked together into a well-ordered 
state, either aligned (ferromagnetic) or anti-aligned (antiferromagnetic), but whose direction 
can rotate. For common forms of magnetic anisotropy the magnetic vector has two or more 
low-energy directions. Several recent investigations have focussed on the possibility that 
the magnetic vector will pass between these directions by quantum tunneling |],^|]. While 
the tunneling rates are predicted theoretically to be exponentially small in the size of the 
magnetic particles, recent experiments ||[4| suggest that quantum tunneling is observable 
in particles with several thousand spins, which are available to present technology. 

In this letter we show that for a wide range of systems, quantum tunneling is completely 
suppressed if the total spin of the magnetic particle is half-integral but is allowed in integral- 
spin particles ||. An important experimental implication of this result is that in ensembles of 
magnetic particles in which the exact number of spins per particle is not precisely controlled 
(the typical case with present technology), half of those particles (those with an odd number 
of electrons) will not exhibit quantum tunneling. Such parity effects are well known in 
atomic physics, but have not been previously noticed for magnetic particles |J . We show in 
several specific examples below that this suppression has a topological origin and arises as a 
destructive interference between different tunneling paths. Thus we find that one quantum 
effect (tunneling) is suppressed by another (interference), leading to "classical" behavior in 
half-integer spin systems! 

We begin with a re-analysis of the tunneling behavior of a small ferromagnetic particle 
with easy-plane-easy-axis anisotropy in the superparamagnetic limit where it behaves like 
a single large spin — referred to as "Model I" in 0. The classical anisotropy energy E has 
the form E(n) = E(6, <fi) = K z o 2 cos 2 6 + K y a 2 sin 2 6 sin 2 (ft, where n is the magnetization 
direction and K z > K y > 0; this energy corresponds to a quantum spin Hamiltonian, 

H = K z a 2 z +K y at (1) 
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where a is the particle's total spin, and Mq = [ibo is its magnetic moment, with \xb the Bohr 
magneton. We are interested in the tunneling of the magnetization direction n between its 
two equivalent low-energy directions at the points 9 = 7r/2, = 0, and 9 = 7r/2, <p = n, 
corresponding to the coherent-state kets |0 > and | vr >, resp. To compute the tunneling 
rate P || we consider the imaginary time transition amplitude expressed as a coherent-state 
path integral for spins |H[H> 
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where f3~ l is the temperature, DQ ~ Y\ T d(j) T d9 T sin 9 T , and where the integral runs over all 
paths (i.e. magnetization directions) connecting one minimum to the other. The Lagrangean 
L occurring in the Euclidean action Se = Jo drL is given by 

L = iacj) — ia<j) cos 9 + E = ia<p + L . (3) 

The first two terms of Eq. (|3]) define the Wess-Zumino term |TIJ which is of crucial impor- 
tance in the following. For any path on the 2-sphere S2, parameterized by 0(r) and 9{r), 
this contribution to the action is equal to ia times the area swept out on S2 between the path 
and the north pole; for closed paths this has exactly the form of the Berry phase flT2]-|T4f . 

The first term of Eq. @ has some special features which require discussion. It is a total 
derivative term, which, when integrated, gives the boundary contribution 

rP 

ia dr<j) = ia((j)((3) - 0(0) + 2mr). (4) 







to Se- Here n is a winding number counting the number of times which the path wraps 
around the north pole. As a total derivative, this term has no effect on the classical (Bloch) 
equations of motion, which can be derived by extremizing just the L piece of the action [0. 
In consequence, this boundary term is commonly ignored. We show here, however, that this 
term is crucial for the quantum properties of the magnetic particle, making the tunneling 
behavior of integral and half-integral spins strikingly different. 

This result is most clearly seen by treating the transition amplitude of Eq. within 
an instanton approximation (i.e., saddle-point evaluation) ||. In the easy-plane-easy-axis 



model, the saddle point paths remain near 6 = ty/2 (cf. [I5H), so they can be characterized by 
their winding in the <fi variable only This model has the important feature that the passage 
from cf) = to <p = 7T can be accomplished either by an instanton or an anti-instanton 
path, i.e., a clockwise or counterclockwise winding over the barrier; see Fig. |l|(a). Then 
the propagator of Eq. (0) is approximated by a sum over paths comprising a sequence of 
instantons and anti-instantons winding over the barrier. The expression for low temperatures 



is 
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Here D is the fluctuation-determinant (without zero- mode) || , I is the number of instantons 
and m the number of anti-instantons in the path, Eq is the zero-point energy in one well, 
exp(— Sq) is the instanton contribution to the path integral, and the constraint on the sum 
reflects the requirement that <j>{0) = ir{mod{2i[)). Note that the action Sq is exactly the 
same for the instanton and anti-instanton, because this action is unchanged if <fi is replaced 



by —cf) [15]. From Eq. (|5|) we can now read off the tunneling rate P (energy level splitting): 



P = 4D\cos(na)\e- s °, e~ s o = (j^^ ) ■ (6) 

Here A = K y /K z . Evidently, the tunneling rate vanishes for half- integral spins because of 
the cos(7rcr)-factor which arises directly from the topological boundary term of Eq. (£|). This 
factor represents an interference between the instanton and anti-instanton contributions to 
tunneling. If the spin a is an integer, then the interference is constructive, and the total 
tunneling rate is of order of the single-instanton rate. But if the spin is half-integral, then 
cos(7rcr) = 0; there is destructive interference between the instanton and anti-instanton, and 
the tunneling rate is zero. Note that this spin-parity effect is of topological origin and thus 
independent of the magnitude of the spin. 

A graphic illustration of the difference between integer and half-integer spin systems is 
afforded by the numerical energy spectra of Hamiltonian ([!]) presented in Fig. |^. Results for 



a = 10|, 10 are shown as a function of K y , for iT 2 = 1. As is increased, a barrier between 
the two easy-axis directions is developed, and tunneling occurs for the integer spin case. The 
tunneling rate, which is proportional to the energy splitting between the two lowest energy 
levels, decreases with increasing K y (see Eq. (^)), eventually vanishing when K y = K z , since 
at this point the x-component of the magnetization commutes with the Hamiltonian and 
so is conserved. In the half-integer spin case, by contrast, it is easy to show from direct 
consideration of Hamiltonian ([[]) that the space of states decomposes into two independent 
subspaces with identical energy spectra. Thus all states are strictly doubly degenerate, and 
there is no tunneling, consistent with the arguments above. 

The result that tunneling is suppressed for half-integer spin has much greater generality 
than is suggested by the above analysis. The suppression can be derived within the coherent- 
state path integral formalism independent of any approximation. Any arbitrary path 9(t), 
4>(t) in Eq. @ (not just the saddle-point path) can always be paired with another path, 
7i — 9(t), — 4>(t), which has the same L in Eq. (|3]), while the winding- number term of Eq. 
(0) is reversed; thus, the destructive interference for half-integral spin occurs term-by-term 
in Eq. (0). Furthermore, this pairing is possible for much more general Hamiltonians than 
Eq. (|I]). One generalization involves adding any additional terms to H which preserve the 
symmetry of rotation around the x-axis (e.g., cr z a y ), including odd-order terms like B ext a x . 
Of course, if such an external-field term is present, then the degeneracy between the two 
wells is broken; nevertheless, the vanishing of the propagator of Eq. (|2|) for half-integer spin 
still implies vanishing hybridization between the two (inequivalent) wells, and an absence of 
tunneling. The other interesting generalization of H involves adding any even-order term 
(e.g., a z , a x a y o 2 z \ but excluding odd-order terms. Then H has time reversal (T) invariance; 
absence of tunneling (i.e., vanishing of the transition amplitude (|5|)) can again be proved 
by a pairing of paths (this time, the paths n(r), — n(/3 — r)). In this case, the vanishing 
of tunnel splitting is directly related to the Kramers degeneracy ||16|| . We can say more 
generally that for anisotropy Hamiltonians which have two equivalent energy minima, the 



absence of tunneling for half-integral spin follows directly from the Kramers theorem |17 



We now analyze a problem raised by a recent experiment on tunneling in small anti- 
ferromagnetic particles 0. The particles were not perfect antiferromagnets, but carried an 
excess magnetic moment, presumed to arise from surface effects, of a few percent of the 



total number of spins. In applying to the experiment the theory developed in Ref. |18| for 
the tunneling of a Neel-ordered antiferromagnetic particle through an anisotropy barrier, 
it was assumed || (as suggested in that the excess spins simply follow adiabatically 

the direction of the Neel vector without affecting the tunneling dynamics. We show now 
that, within a simple extension of the existing models, this expectation is false: As above, a 
half-integer excess spin of any size quenches tunneling; an integer excess spin also influences 
the tunneling rate, though more modestly. 

Our model is a ID antiferromagnetic ring with N spins tr^> (spin magnitude a) and 
periodic boundary conditions; all the spins cr are coupled to a central excess spin s with 
coupling constant J^> = (— l)- 7 J c . Thus the central spin will prefer to remain aligned with 
the Neel vector of the spins on the ring (see Fig. |TJ(b) ) . The Hamiltonian of the system may 
be written 

N 

H = ^[J<T W -<7^ +1 ) + H(tT®) + J®<T®*]. (7) 

3=1 

Here N is even; odd N would produce a different problem involving frustration and non-zero 
total spin on the ring. H is the anisotropy Hamiltonian ([TJ), and leads to two degenerate 
states, the one shown in Fig. |l](b), and the one where all the spin directions are reversed. 
The problem is to compute the tunneling rate between these two states. 

This model may be considerably simplified if we restrict our consideration to the con- 
tinuum and semi-classical (large a) limit. For this we use again the coherent-state path 
integral representation, and, applying standard manipulations we then arrive at a 

generalized non-linear cx-model in terms of a Neel unit- vector 1, and a central spin unit vector 



n s |TjJ. We take 1 to be uniform around the ring, obtaining the following effective Euclidean 
action: 



S e= I drAief + tf sin 2 OD+NEiOi^i) 

JO OUr 



6 



+aNJ c l-n s + is<p s (l- cos6 s )], (8) 

where 1 and n are expressed in polar coordinates, and the transverse susceptibility x± is 
related to the parameters of Eq. (0) by x± — N/a^ / J. 

As before, we find that if the excess spin s of the antiferromagnet is half-integer, then the 
tunneling rate is exactly zero. The proof follows in the same way: we consider an arbitrary 
path {9i(t), <f>i(r), s (t), s (t)} and its partner with 0i s — > n — 9[ s , <p^ s — > — <p^ s . Then in the 
path integral expression analogous to Eq. fl2|), and with 1 and h s having the same boundary 
conditions, these two paths have opposite winding- number contributions ±zs7r(l + 2n) (cf. 
Eq. (|^)), and the same values of L . Thus, a factor cos(7rs) appears exactly as in Eq. 
(||), implying complete destructive interference and a vanishing of the tunneling rate if the 
central spin s is half-integral. 

Also as in the ferromagnetic case, this vanishing can be seen to be related to the Kramers 
degeneracy. Again, the model Eq. (0), and therefore Eq. (f|), has time reversal invariance 
(since all terms contain an even number of spin operators). Thus, the ground state is a 
Kramers doublet so long as the total spin of the model is half-integral, which, since N is 
even, requires that s is half-integral. Again, suppression of tunneling is related to the absence 
of a tunnel splitting in the ground state. However, we caution as we did above that there 
is not a one-to-one correspondence between the Kramers theorem and absence of tunneling 
for antiferromagnetic models. 

Finally, we consider the question of how strongly a nonzero but integer s modifies the 
tunneling dynamics of Eq. (|8|). For s = this model can be solved in the instanton 
approximation; the saddle point solution happens to be identical to a different (uniaxial- 



anisotropy) model considered in JT8|. The saddle point path has 9\ = n/2 everywhere, while 



4>i passes from to it. As in [IS] , the tunneling rate is given by: P ~ uje~^ 2x±ky ^ B , with 
cu = ^B(8k y /x±)K and with the notation k y ^ z = a 2 NK y ^ z for the anisotropy constants of 
the whole ring. 

A full solution of Eq. fl8|) for s ^ seems difficult. We can make progress in the adiabatic 
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approximation , in which the spin n s simply follows the instantaneous direction of the 

Neel vector 1 at every point along the path. In Eq. this simply involves removing the 
1 ■ n s term (since it is just a constant) and setting 9 S = 9i, <p s = (pi ptifl . 

We can now find an approximate solution for the tunneling rate if we assume that k y <C k z 
and that therefore 9i does not fluctuate very far away from tt/2. If we write 9[ = 7i/2 + •§ 
and expand the Lagrangean to second order in we obtain 

-{$ 2 + <p 2 ) + $ 2 {k z -k y sm 2 <f)- " 



+k y sin 2 4> + is<p(l + $). (9) 

We find |ZT] that for k y <C k z we have ft 2 <C <p 2 <C 8fi 2 B k z /x±- Keeping then only the leading 
terms in L we find for the extremal $-path, d ~ —is<j)/2k z , obtaining an effective Lagrangean 
for the variable, L w (xj_/8/i^ + s 2 /4k z )4> 2 + k y sm 2 (ft + is0. This Lagrangean is now 
identical to the one for a pure antiferromagnet (see HT8| ) with s = 0, but with a topological 
boundary term and a modified value of the transverse susceptibility (the "moment of inertia 
of the rigid rotor" |IT| ): x± = X± + 2s 2 /i^//c 2 . Thus the tunneling rate is finally given by 

P ~ |cos(tts) Icue-^xf^y/^ . (10) 

Again, for half-integral excess spin s tunneling is suppressed. For integral s we can define 
a crossover excess spin s c for which the contribution of the excess spin to the moment of 



inertia becomes comparable to that of the antiferromagnetic ring: s c = \Jx±k z /2/ fig. If, 
e.g., k z m lQky, then s c is approximately 1.6 times the WKB exponent of Eq. (|T0D for s = 0. 



As Ref. |18j has pointed out, for practical reasons this WKB exponent cannot be much larger 
than about 25 if tunneling is to be observed in a small particle. Thus, the magnitude of 
the excess spin is quite restricted (s < 1.6 x 25) if pure antiferromagnetic dynamics is to be 
observed. 

In summary, we have demonstrated that tunneling in a wide class of magnetic particles is 
strongly parity dependent, being completely suppressed for half- integral spins. Preliminary 
results indicate that similar effects are to be expected in the tunneling of domain walls [H. 



We expect that these phenomena can still exist in the presence of moderate dissipation [11 
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FIGURES 

FIG. 1. (a) Anisotropy energy vs. at (9 = ir/2, showing the path for the instanton (T), the 
anti-instanton (AT), and a more general path V containing one anti-instanton and four instantons. 
(b) Antiferromagnetic ring coupled to an excess spin (Eq. (|7|)). 

FIG. 2. Low-lying eigenenergies as a function of K y for K z = 1 for the model of Eq. (|l|) for 
spin-10| and spin-10. All levels are doubly-degenerate for a = 10^. 
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